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Abstract. We review the techniques used to renormalize quantum field 
theories at several loop orders. This includes the techniques to sys- 
^ ■ tematically extract the infinities in a Feynman integral and the imple- 

£ — ' mentation of the algorithm within computer algebra. To illustrate the 

CO , method we discuss the renormalization of (f> 4 theory and QCD including 

the application of the critical point large N technique as a check on the 
anomalous dimensions. The renormalization of non-local operators in 
QCD is also discussed including the derivation of the two loop correction 
to the Gribov mass gap equation in the Landau gauge. 
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1 Introduction 



Our theoretical understanding of quantum phenomena in particle physics and 
condensed matter is guided by the underlying quantum field theory. This is for- 
mulated in terms of a Lagrangian of fields which have symmetries motivated by 
experimental observation. For instance, the electromagnetic field governing light 
leads us to the gauge principle and gauge field theories, and in particular quantum 
electrodynamics. The properties of this field theory have been generalized to the 
current theory describing all elementary particles which is the standard model. In 
order to determine predictions from these quantum field theories requires one to 
develop the loop expansion of perturbation theory where the interactions of the 
fields are represented by Feynman diagrams. One of the main properties of such 
diagrams is that when one evaluates them they are divergent. However, as has 
been established for a long time there is a systematic and mathematical procedure 
for handling the resulting infinities which goes under the general title of renormal- 
ization. For example, a comprehensive survey of the area is provided in [J. The 
topic is by nature a technical one and if one wishes to extract meaningful and ac- 
curate predictions from a quantum field theory one needs to be able to perform 
the renormalization at a large number of loops. In this article we review some of 
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the practicalities of achieving this in several field theories which are of main inter- 
est. These are scalar <p A theory which is relevant to phase transitions in various 
spacetime dimensions and quantum chromodynamics (QCD) which is the quantum 
field theory describing the quarks and gluons associated with the strong nuclear 
force. Our aim is to address the basics in the first part of the article where we 
discuss various techniques for evaluating the complicated and divergent Feynman 
integrals which arise in renormalizable field theories, and the way they are handled 
in practice by using computer algebra and symbolic manipulation packages. Since 
the renormalization is a complicated procedure we also indicate the rudiments of 
the important task of how one actually verifies that a calculation is in fact cor- 
rect. The second part of the article discusses recent problems of interest which 
apply the techniques discussed in the first part and illustrate their application from 
a practical point of view. These will primarily be in the context of QCD where, 
for instance, its three loop renormalization is discussed in a particular non-linear 
gauge. We will also consider the problem of renormalizing operators which have a 
particular degree of non-locality. 

The article is organized as follows. We outline the main issues concerning 
renormalization in section 2 before discussing various calculational techniques in 
section 3. Section 4 surveys the main checks one has on the derivation of the 
renormalization group functions which lead on to more specific checks for QCD in 
the large Nf limit in section 5. We discuss more recent aspects of renormalization 
of QCD in section 6, before examining the issue of how one can treat a particular 
class of non-local operators in section 7. Finally, various concluding remarks are 
provided in section 8. 

2 Statement of problem 

We commence by summarizing from a general point of view the main issues 
underlying the renormalization of a quantum field theory and the main terminology 
of the subject used in this article. At the outset it is worth doing this with several 
basic field theories in mind but initially we will focus on the scalar quartic interac- 
tion, cf> 4 , in four spacetime dimensions. Later we will consider theories with gauge 
symmetries and in particular QCD. However, the general remarks and comments 
will apply equally to all field theories which are renormalizable in four or other 
dimensions. For <jr theory the basic Lagrangian is 

L = \ (dcbof + |^ (2.1) 

where the subscript denotes that the field (f> and coupling constant g are bare 
quantities. This observation is essentially the foundation of the problem of the 
need for renormalization. If one were to naively compute quantum corrections to 
any Green's function with (|2.1|l as the starting point, then the resulting Feynman 
diagrams would be infinite in four dimensions. This is not surprising as one is using 
a (local) Lagrangian involving quantum fields defined at the same spacetime point. 
In other words the variables, parameters or fields of the initial Lagrangian can be 
regarded as being insufficient or not the correct variables in which to define the 
problem. To circumvent this, one defines a new set of (renormalized) variables by 
scaling the original variables with a (multiplicative) factor, known as the renormal- 
ization constant. For l|2.1|) we define these formally as <j> — and g Q = Z g g 
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which leads to the Lagrangian in terms of renormalized variables 

L = {d<t>f + { x Z\Z g ^ . (2.2) 

This establishes the framework for renormalization. If one subsequently computes 
the previously infinite Green's functions with these, as yet undetermined renormal- 
ization constants, then by choosing their value appropriately the infinities can be 
absorbed into the renormalization constants, Z$ and Z g . As the quantum field 
theory is renormalizable in four spacetime dimensions then the choice one makes 
for say the divergent 2-point and 4-point Green's functions of c/> 4 theory, means 
that all other n-point functions are finite. Non-renormalizable theories require ad- 
ditional operators in the Lagrangian over and above the original ones to retain 
finiteness. By contrast a superrenormalizable theory either does not require all the 
available renormalization constants or these actually have a finite number of terms 
when expanded in a coupling constant expansion. Whilst this is the general proce- 
dure, one ordinarily fixes the renormalization constants order by order in a coupling 
constant or perturbative expansion. Though one can still follow this prescription 
non-perturbatively such as in a lattice regularization. 

Whilst this is the overall essence of renormalization, there are clearly several 
technical issues to be addressed which we now briefly discuss. The choice in redefin- 
ing the variables may seem ad hoc and has little connection with the real physical 
world. Indeed one may not only absorb the infinities of the Green's functions but 
also an arbitrary finite part. How one practically removes the divergences is known 
as the renormalization scheme. Although there are a large number of such schemes 
the most widely used is the modified minimal subtraction scheme denoted by MS. 
To ensure that the same physics emerges independently of how the infinities are ab- 
sorbed, the information contained in the renormalization constants are encoded in 
renormalization group functions such as 7^(3) and (3{g) which appear in the renor- 
malization group equation. These functions determine properties of the quantum 
theory. For instance, in field theories underlying condensed matter problems the 
critical points relating to the phase transitions are given by the non-trivial values 
of the coupling constant, g c , where (3(g c ) = 0. In particle physics applications, for 
example, the solution of the differential equation defined by the /3-function deter- 
mines how the running coupling constant depends on the renormalization scale, 
which is defined later. For instance, in QCD since the /3-function is negative for 
small values of the coupling constant then the theory is asymptotically free, |2j. In 
other words the (confined) quarks and gluons behave at very large energies as if 
they were free particles. 

In alluding to infinities in general terms so far, for practical calculations one 
must have a mathematical way of handling them and determining their nature. 
Moreover, for theories with a classical symmetry which is assumed to be preserved 
in the quantum theory, the procedure for quantifying or regularizing the infinities 
must respect the symmetries of the theory. In this context we briefly mention the 
algebraic renormalization technology developed primarily to determine the general 
consequences for the renormalization of gauge theories but also for theories with su- 
persymmetry. For example, see 0- In essence algebraic renormalization determines 
the form of the renormalization constants consistent with the Slavnov- Taylor iden- 
tities. Two of the most widely used regularization procedures are lattice regular- 
ization and dimensional regularization. In the former one replaces or approximates 
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the continuum of spacetime by a discrete set of points regularly spaced a distance 
a apart 1 . The presence of a is a feature of any regularization which is that every 
regularization introduces an arbitrary scale. For lattice regularization one wishes 
to take a as small as possible in order to be as close as possible to the continuum. 
Though in practical (and financial) terms this is very costly from the point of view 
of computer running time. However, the point is that whatever physical quantity 
one is computing the result has to be independent of the arbitrary regularization 
scale. 

The same feature is true in the other regularization we consider. By contrast, 
one retains a continuum spacetime but changes the dimension away from four by 
an infinitesimal amount by defining the spacetime dimension as d = 4 — 2e where 
e <C 1 and e > 0. In this regularization the infinities in Green's functions manifest 
themselves as poles in e and after the poles are removed by some criterion in some 
scheme, then the e — > limit can be taken in a non-singular way. However, in 
changing the spacetime dimension an arbitrary scale is introduced by requiring 
that the coupling constant remains dimcnsionless in rf-dimensions. We now take 

g a = fj, 2e Z g g (2.3) 

where [i is the arbitrary renormalization scale. Since the Green's function after 
renormalization, T^ n \jjL,p,g, . . .), will necessarily depend on fj, in the combination 
p 2 //x 2 where p is a momentum, we can quantify the essence of the renormalization 
group as 

M frW( M ,p, 9 ,...) = (2.4) 

so that results are independent of the arbitrary scale /x. 

One point of clarification needs to be made in the context of defining (dimen- 
sional) regularization which is that infinities fall into two classes. Those which 
are ultraviolet arising from divergences in Feynman integrals at large momenta and 
those which are infrared coming from divergences in integrals at low loop momenta. 
In massless theories dimensional regularization regularizes both types of infinities 
so that it is never clear where the poles in e originate from. However, infrared infini- 
ties can be regularized by a non-zero mass which clearly acts as a low momentum 
cutoff. This mass can either be present in the original Lagrangian or introduced by 
hand. In the latter case when the theory is renormalized it can be smoothly set to 
zero. However, in gauge and supersymmetric theories such extra masses can break 
the symmetries and if such regularizing masses are introduced, care is required to 
preserve the symmetry. 

We close this section by briefly summarizing several renormalization schemes 
which arc in common use. We have alluded to MS already and we note that the 
original minimal subtraction scheme, MS, on which it is based requires the removal 
of the poles in e only. (Though it is worth noting that MS and MS are not solely tied 
to dimensional regularization. One can minimally subtract divergences in lattice 
regularization.) The MS scheme is a modification of MS in that a common finite 
term, Ane~^ where 7 is the Euler-Mascheroni constant, is removed in addition to 
the poles as it was observed that in practical calculations that the convergence 
of perturbative expansions was improved, 0]. These two schemes fall into the 
class of mass independent renormalization schemes. By contrast mass dependent 



^his is known as a regular or square lattice which is the most widely used. However, lattice 
regularization is not restricted to square lattices. One can define triangular or random lattices. 
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schemes are those where, aside from the poles in e being removed certain finite 
parts are also removed which involve masses of the fields introduced by putting the 
external particle masses on-shell or setting a Green's function to a particular mass 
independent number at some value of the renormalization scale. Examples of such 
mass dependent schemes are MOM, [5], on-shell and the modified regularization 
(RI'), 0, schemes. In lattice regularization RI' is one of the main choices as it is 
constructed in such a way as to minimize computer calculation time. In a similar 
way to requiring that physical quantities are independent of the renormalization 
scale introduced through the regularization, they also have to be independent of 
the renormalization scheme. However, it is possible to convert between different 
renormalization schemes. Indeed the MS scheme is used as the reference in this 
context primarily because one can compute much further in the coupling constant 
expansion than in a mass dependent renormalization scheme. Therefore, on the 
assumption that with more terms the series is more accurate, it should be closer to 
the real physical value of the quantity being computed. So, for example, in lattice 
computations which are non-perturbative but computed in the RI' scheme, one has 
to convert the results to the MS scheme and match the large momentum behaviour 
to the perturbative MS result at whichever perturbative order is available. Thus 
the larger the number of terms in the coupling constant expansion that are available 
then this will reduce the error in the final numerical estimate. For certain problems 
this has been achieved at three and four loops, [7||S]. 



Having discussed the method of renormalization in general terms, we now de- 
tail some of the techniques used to extract regularized infinities from Feynman 
integrals. We will concentrate on the application of dimensional regularization to 
massless and massive Feynman diagrams. Throughout we will regard the integrals 
as expanded near four dimensions with d = 4 — 2e. First, we recall a simple but 
powerful technique to treat one loop diagrams which is based on a simple inte- 
gral and attributed to Feynman. Within a Feynman graph one has products of 
propagators and the basic idea is to write these as an integral using 



where notionally a = k 2 — m 2 and b — (k — p) 2 — m\ with k regarded as an internal 
momentum and p as external. With the two momenta within one factor, one can 
use Lorentz symmetry before performing the fc-integration. This leaves a function 
of the Feynman parameter, x, which at one loop can be related to known functions. 
For example, when the masses mi and m,2 are equal the integral can be written 
in terms of a hypergeometric function. Whilst this is a powerful approach at one 
loop, it ceases to be practical at higher loops especially when there are more than 
one mass and external momenta. In practice one chooses a method of calculating 
the Feynman integral which is tailored to the overall (renormalization) problem of 
interest. 

We will focus on one of these and consider the situation where the field theory 
involves particles of one non-zero mass and is renormalizable. Then to renormalize 
the divergent Green's functions and to extract the infinities, one expands the set 
of (scalar) Feynman integrals in powers of p 2 where p is the external momentum. 
This method is known as the vacuum bubble expansion. In regarding all integrals 
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as scalars, we have assumed that one has first decomposed Lorentz tensor integrals 
into their scalar pieces. To focus on the particulars of the vacuum bubble expansion 
we consider the simple example of a self energy type integral in a 2-point function 
which gives 



1 f 1 f p 



,-.2 



[k 2 — m 2 ][(k — p) 2 — m 2 ] J k [k 2 — m 2 ] 2 J k [k 2 

4 



In the final term the (kp) 2 numerator factor has been simplified by using Lorentz 
symmetry. The expansion truncates due to the renormalizability condition which 
states that since the theory is renormalizable the 0((p 2 ) 2 ) terms are finite since 
otherwise one would require a fourth derivative 2-point term which is forbidden by 
renormalizability. The resulting integrals on the right hand side of l|3.2|l are simple 
vacuum loops which can be evaluated. More appropriately the method can be easily 
extended to higher loops. For instance, in four dimensions three loop single mass 
scale vacuum bubbles are known to their finite part and at two loops, three mass 
scale vacuum graphs are also known to their finite part with respect to e. See, for 
example, [jj]. One of the benefits of this algorithm is that it can be implemented 
in computer algebra in an automatic way. Also the method is applicable to higher 
leg Green's functions if one wants to renormalize them too. However, in certain 
problems where the finite part is required exactly in, say, scattering problems, the 
more recent technique of is appropriate. There each propagator is represented 
by a Mellin-Barnes integral and contour integration used to evaluate the two loop 
4-point functions, for example, to their finite parts. Though currently results have 
yet to be derived for the physically relevant case of all possible massive propagator 
combinations. 

Next, we note that one distinct advantage of the massive vacuum bubble expan- 
sion is that there are no infrared problems as the inherent mass provides a natural 
infrared regulator. For theories where there are massless fields it may seem that 
the vacuum bubble approach is inapplicable. However, one can manufacture 
a fictitious mass p, which acts as an infrared regulator via the identity 

1 = 1 [2kp-p 2 -fl 2 } 
{k-p) 2 [k 2 -p?] (k - p) 2 [k 2 - Ji 2 \ ' 1 ' ' 

Within a Feynman diagram this identity can be used repeatedly with the truncation 
criterion based on Weinberg's theorem, for the overall finiteness of a Feynman 
integral. Like the completely massive case this algorithm can be implemented 
automatically in computer algebra. Though one ought to be aware of the potential 
problem of breaking an inherent symmetry of the theory when a non-zero /2 is 
introduced. For non-abelian gauge theories this is discussed in 

Another equally powerful approach for massless field theories, such as QCD, is 
the use of integration by parts based on the identity, in d-dimensions, |13j . 

= J 'w?w^ I{ »- k >i < 3 ' 4) 

where /(p, fc, ....) is the integrand derived from the propagators and vertices. Al- 
though such identities can be used for massive integrals, the power lies in the 
observation that the differentiation can introduce numerator propagator factors to 
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simplify the momentum topology. By momentum topology we mean that diagram 
which represents all the scalar propagators. This is not necessarily the same as 
the actual Feynman diagram topology itself as the cancellation of a denominator 
factor means that that line would be omitted in the momentum topology. This 
cancellation is one of the principles which underlies the Mincer algorithm, |14) . 
This is a package for evaluating massless three loop 2-point functions in dimen- 
sional regularization to their finite parts. Moreover, it has been encoded in the 
symbolic manipulation language FORM, llfij . Essentially at three loops there 
are fourteen basic integration topologies each with its own integration by parts 
routine. Though these momentum topologies are not all independent. Whilst pri- 
marily used for 2-point functions, Mincer can be applied to 3-point functions when 
one of the external momenta is set to zero or nullified. However, in this case one 
must ensure that the nullification does not introduce spurious infrared singularities 
such as J k ^2)2 which result in poles in e which being infrared in nature, cannot 
be distinguished from ultraviolet ones. Despite this technicality there are methods 
for handliing nullification known as infrared rearrangement, |17| . Though this has 
not been implemented automatically in computer algebra. The main application of 
Mincer is to the renormalization of four dimensional gauge theories at three loops 
and in particular QCD. 

Given this emphasis on automatic computer algebra we make some specific 
remarks. The need for such machinery can be gauged from the fact that when 
one increases the loop order of a calculation, the number of Feynman diagrams 
increases almost exponentially. As examples we note that the recent full three 
loop renormalization of QCD in the maximal abelian gauge required the evaluation 
of 37322 diagrams, Also the four loop QCD /3-function of in a linear 

covariant gauge required of the order of 50000 diagrams. Clearly computers are 
necessary to not only implement the computational algorithm, such as Mincer or 
the vacuum bubble method, but also to handle the sum of the individual results. 
Packages such as Mincer, which have been optimized, are essential to having as 
short a computation time as possible. Indeed a publicly available four loop Min- 
cer would be equally as useful. For such automatic computations the Feynman 
diagrams themselves need to be generated electonically and Qgraf, (SHj, has been 
developed specifically for this purpose. It has various output formats which can be 
readily converted to the notation used by Mincer before applying the algorithm 
itself. The implementation of the renormalization procedure can also be performed 
automatically without the need for the traditional method of subtractions. This 
method determines the absolute divergence of a diagram by subtracting all sub- 
graph divergences. Instead in the automatic approach, |21| . the Green's function is 
computed as a function of the bare parameters. Then the counterterms are intro- 
duced naturally (and equivalently to the subtraction method) by rescaling by the 
perturbatively expanded renormalization constants g = Z g g and so on. Once the 
counterterms have been implemented at a particular loop order, the divergence re- 
maining is then that associated with the renormalization constant of that particular 
Green's function. 

Finally, we briefly comment on other recent techniques of evaluating Feynman 
diagrams. One which is also appropriate to n-point functions at two loops for non- 
zero external momenta is the differential equation method of [221- The basic idea 
is to derive a complete set of differential equations at a particular loop order for 



8 



J. A. Gracey 



relevant momentum topologies. These are then solved with a basic set of master 
integrals as the boundary conditions which have to be evaluated by direct meth- 
ods. For completely massless integrals in problems where conformal symmetry is 
present, such as at a fixed point, the method of uniqueness, [2HII2], is also power- 
ful. Essentially when the sum of powers of three momenta in a loop integral satisfy 
a particular condition depending on the spacetime dimension, then the integral can 
simply be replaced by products of related propagators and Eulcr T-functions of 
the original propagator exponents. For large Nf methods, which will be discussed 
later, this has proved to be extremely powerful in computing information on the 
renormalization group functions beyond the leading large Nf order and to all orders 
in perturbation theory. 



Performing a renormalization even of a simple quantum field theory to several 
loop orders can lead to the computation of a large number of Feynman diagrams 
which are determined by one or other of the methods previously discussed. However, 
one natural question arises in the ultimate derivation of the renormalization group 
functions and that is whether the results are correct. For major loop calculations in 
four dimensional gauge theories to check by repeating the work using independent 
computer algebra programmes will not only be time consuming but make a large 
demand in both human and computer resources, which could be used for other 
problems. The exception to this is the situation where one is extending existing 
anomalous dimensions and /3-functions to the next loop order. This requires the 
counterterms which depend on the finite part of the Green's functions which were 
renormalized in the previous loop calculation but which are not ordinarily computed 
to construct the established renormalization constants. Therefore, one invariably 
has to reconstruct the previous calculation prior to tackling the new loop order. 
Aside from this check we now address one standard way of assessing whether the 
calculation one has performed is correct. By this we mean that it satisfies a set of 
internal consistency checks at the very least and others tailored to the problem in 
hand. First, the renormalization group provides a clear insight into the structure 
of the renormalization constants. To illustrate this, if we work from the form of the 
MS (3- function and anomalous dimension 7^(3) given by 



in d = 4 — 2e dimensions where 3 is a generic coupling constant then the corre- 
sponding renormalization constants must be of the form 



4 Checks 



0(9) 
14(9) 



(d - 4) 5 + Ag 2 + Bg 3 + Cg 4 + 0{g b ) 
ag + bg 2 + eg 3 + 0{g 4 ) 



(4.1) 



Z 
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(d-4) 
A 3 



TAB 



C 



9 3 + 0(g 4 ) 



+ 



(d-4) 3 



+ 



6(d-4) 2 



3(d-4) 




[ (2a A 2 - 3a 2 A + a 3 ) (3a6 - 2aB - 2b A) 
+ [ 6(d-4) 3 + 6(d-4) 2 
+ 0(g 4 ) 



3(d-4) 



c 



(4.2) 
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to three loops. Clearly the residues of the simple poles are in a one to one cor- 
respondence with the coefficients of the renormalization group functions since the 
MS scheme is used. However, the residues of the other poles depend only on the 
structure of the previous loop calculations. Therefore, in a new loop calculation 
these poles are already predetermined and must emerge from the new one loop cal- 
culation. This is important in automatic calculations since one determines the full 
renormalization constant without using the subtraction method. Another internal 
check is that provided by symmetries of the original theory. Ignoring the techni- 
calities produced by anomalies, in a gauge theory certain renormalization group 
functions are independent of the choice of gauge in mass independent renormal- 
ization schemes. Therefore, working in arbitrary covariant gauges means that the 
gauge parameter must be absent in the final result. By contrast, when renormaliza- 
tion group functions which depend on the gauge parameter, a, are calculated, the 
residues of the triple and double poles of the renormalization constants will depend 
on a and also be constrained by the conditions above. Briefly, in supersymmetric 
theories similar conditions emerge. For example, when supersymmetry is unbroken, 
the anomalous dimensions of the fields in the same supermultiplet have to be equal 
when the component field version of the Lagrangian is considered. This is a strin- 
gent check not only on the renormalization procedure but also on the regularization 
which must preserve the supersymmetry. Whilst this ensures the renormalization 
constants are checked to an extent, the residues of the simple poles are not. How- 
ever, there are partial checks available in some theories from the large Nf methods 
developed for scalar field theories, [241 I25ll2l)| . 

This can be illustrated in the case of 0(N) 4> 4 theory where the scalar field is 
a vector in O(N). In MS in d-dimensions the /^-function has the form 



/%) = (d-4)| + [N + %] 9 - - [3N + U]L. 

o 4 

+ [337V 2 + 922/V + 2960 + 96(5/V + 22)C(3)] j— + 0{g 5 ) (4.3) 



where C( n ) is the Riemann zeta function. In e?-dimensions there is a non-trivial 
fixed point of the /3-function, g c , defined by f3(g c ) = known as the Wilson-Fisher 
fixed point. Clearly g c will be a function of d and N. Though near four dimensions 
with d = 4 — 2e we have g c — g c (e, N), which can be expanded in powers of 1/N 
as N — > oo giving 



Evaluating the renormalization group invariant universal critical exponent (3'(g c ) 
in the same limit gives 



which encodes the information of the original /3-function, albeit in a different way. 
However, if one can compute f3'(g c ) in d-dimensions order by order in the 1/N 
expansion then the coefficients in the polynomials of N in 14.3(1 can be read off 
directly. In [21 12H 1251 I2Z| the d-dimensional large N technique was developed for 
the O(N) non-linear a model in d = 2 + e dimensions and the first three terms in 
the 1/N series for the exponents rj, v and u> have been determined in (^-dimensions. 



.2 
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As this model is in the same universality class as O(N) 4> 4 theory at the Wilson- 
Fisher fixed point then one can partially check off the known coefficients in the 
corresponding renormalization group functions. For instance, writing the leading 
large N wave function anomalous dimension as 



d r N 



• e r )N r - 2 g r+1 



then with U 6 ,2 = T,n> m >o 
at ten loops and 0(1/ N 3 ) that 



7 ( 5 ) = ]T(c r iV 2 

as the (3, 4) torus knot number, 



r=l 

c-ir 



(4.6) 
J, one finds 



e 9 = [1560674304C(10) - 12534896640C(9) + 11070010560C(8) 

+ 1732018176C(7)C(3) + 581961984C(7) - 3411394560C(6)C(3) 

- 2684240640C(6) + 209534976C 2 (5) - 1567752192C(5)C(4) 
+ 1754664960C(5)C(3) - 975533568C(5) - 9289728((4)C 2 (3) 
+ 1310201856C(4)C(3) + 1636615872C(4) - 137158656C 3 (3) 

- 1708996608C 2 (3) + 294403968C(3) 

- 89800704f/ 62 - 34 1 35 4 33] /1950396973056 . (4.7) 

This requires not only knowledge of r\ at 0(1/ N 3 ) but also w at 0(1/ N 2 ) as this 
encodes the value of g c required for rj = "f(g c )- By way of illustration as to the form 
of such large N exponents we note that with 



LO = fl — 2 



E 



iV 7 



then 



and, [27], 



UJ 2 



uji = (2/i - l) 2 r/i 
4( M 2 - 5 M + 5)(2n - 3) 2 (/x - !)//(#(// 



(4.8) 



(4.9) 



* 2 (//) 



(/i-2)3( M -3) 



16/i(2/i - 3) 2 

( M - 2)3(^-3)^! 

3(4// - 48// + 241 (i 3 



549// + 566/i - 216) (/z - l)/x 2 0(/i) 



2(^-2)3(^-3) 
[16//° - 240// + 1608// - 6316// + 15861// 
25804// + 26111/ - 14508// + 2756// 
672 M -144)]/[( M -2) 4 ( M -3) 2 ]*(/i) 
[144// 4 - 2816// 3 + 24792/Z 2 - 130032/i 11 + 452961//° 
1105060// + 1936168/t 8 - 2447910// + 2194071// 
1320318// + 460364// - 43444// - 26280// 

8208/; - 864]/[2(2/t - 3)(/x - l)(/t - 2) 5 (/i - 3) 2 //] 1 yy 2 (4.10) 



where r/i = - 4r(2/i - 2)/[r(2 - //r(/i - l)r(/x - 2)r(/i + 1)] and we have set 
d = /t/2. The expression for 102 contains derivatives of the T-function which are 
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denoted by 

#(/i) = ^(2 M -3) + ^(3-/x) - - </>(!) 

e(/x) = v'(m-i) - ^'(1) 

$( M ) - ^'(2M-3) - " ^'(M-l) + (4-11) 

Whilst 0(N) cj) 4 theory has been examined to several orders in the same 

technique has been developed for QCD in [2131211 • There the expansion is in terms 
of the number of quark flavours, Nf, and not N c which is the number of colours. The 
latter 1/N C expansion deals with the structure of QCD from a completely different 
point of view. 

5 Large N f QCD 

Having outlined the critical point approach for O(N) 4 theory, we now con- 
centrate on how one develops the same formalism for QCD in the large Nf limit. 
Examining the d-dimcnsional MS j3- function for QCD, one finds that there is an 
equivalent Wilson-Fisher fixed point which can be accessed in powers of 1/-/V/ as 
a function of e. However, to compute critical exponents in QCD in d-dimensions 
and relate them to the associated anomalous dimensions in perturbation theory at 
whatever order they are available, one does not compute with the QCD Lagrangian 
itself. Instead one exploits the properties of the d-dimensional fixed point by realis- 
ing that at this fixed point QCD is in the same universality class as the non-abelian 
Thirring model (NATM), . This is a four-fermi theory which is renormalizable 
in two dimensions and plays the role of the non-linear a model in the previous 4> A 
theory example. The Lagrangian of the NATM is 

L NATM = iQUfrpI + ^^tfj^yyj)* (51) 

where A is the coupling constant of the NATM and is dimensionless in two di- 
mensions. Introducing an auxiliary spin-1 field the Lagrangian can be written as 

A a2 

L NATM = fyilfyil + A a ^ lI Tj a j^llj lJ - . (5.2) 

Clearly there is no field strength term in the non-abelian Thirring model and from 
a critical point of view this is due to the fact that at this infrared stable fixed point 
that operator is irrelevant. The critical behaviour is driven by the common quark 
gluon vertex when one compares with the QCD Lagrangian 

L QCD = _ \ Q a naiiv _ J_( 9 M^«)2 _ c a d"D^C a + i$ U fllp U (5.3) 

4 2a 

where in this section we use a as the covariant gauge parameter with a = cor- 
responding to the Landau gauge. Within the large Nf critical point formalism the 
triple and quartic gluon vertices are embedded inside quark loops with respectively 
three or four external auxiliary spin-1 fields, |31| . 

As an example of the large Nf critical point formalism, the 0(1/Nf) correction 
to the /3-function can be computed by considering the critical point structure of 
the gluon 2-point function with the insertion of the operator G a ^G a ^ v : 29 . The 
anomalous dimension of this operator will be related to the critical exponent of the 
associated coupling constant which is clearly the QCD /3-function when expanded in 
powers of e near four dimensions since u> = — f3'(g c ). At criticality the propagators 
of the fields of the non-abelian Thirring model or QCD Lagrangian obey asymptotic 
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Figure 1 Feynman diagrams for 0(1/ Nj) correction to w. 



scaling forms where the canonical exponents are derived from the dimensionality of 
each term in the d-dimensional Lagrangian. Thus we have, in the Landau gauge, 

US EDI, 



k 2 



(5.4) 



where A and B are momentum independent amplitudes and 

a = M- l + |»7,/3 = l- »?-X- (5-5) 

The exponents r\ and x are respectively the quark anomalous dimension and the 
quark gluon vertex anomalous dimension at criticality. In large Nf there are four 
diagrams contributing to the gluon 2-point function with the operator inserted 
which are illustrated in Figure 1. The three loop diagram is included since it is of 
the same order in large Nf as the two loop ones because in the large Nf counting 
one regards the amplitudes as A = 0(1) and B = 0(1 /Nf). These diagrams are 
infinite but can be regularized by shifting /3 to /3 — A where A is the regularizing 
parameter, E31 HH1 HHI- It plays a role akin to e in dimensional regularization 
and it is important to appreciate that within the critical point large Nf method 
e is not the regularization. The calculations are performed in fixed d-dimensions. 
After subtraction of the poles in A, using the renormalization procedure of \',V2\ . 
the scaling behaviour of the remaining finite Green's function determines u> and we 
have, ESI, 



(fx-2) - [(2(i-3)(fi-3)C F 

(4/x 4 - 18^ 3 + 44^ 2 ~ 45^ + U)C A 



(5.6) 



T(R)N } 



4(2/x-l)( M -l) 

where r) = (2/x - 1)(/j, - 2)r(2^)/[4r 2 (^)r(^ + 1)T(2 - fj,)]. The factors deriving 
from the group theory are defined by 

T a T a = C F I , f acd f bcd = C A S ab , Tr (T a T b ) = T F 5 ab . (5.7) 
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Consequently one can verify that the known leading 1/Nf coefficients agree with 
the MS four loop /3-function of |19j . Moreover, the unknown higher loop coefficients 



at 0(1/Nf) are encoded in 
of the QCD /3-function as 



. For instance, if we write the leading large Nj part 



P(g) = fasf + ^a n+1 [T F N f ] n g 



nn+2 



1 



(5.8) 



then at five loops we have 

a 5 = [(288C(3) + 214)C f + (480C(3)- 229)Ca]/31104. (5.9) 

For O(N) <t> A theory, it was possible with the critical point method to go to 
several orders beyond the leading 1/N order. For QCD the same is possible, in 
principle, but one has the added complication of the Lorentz structure associated 
with fermions and the fact that the structure of the quark gluon vertex is such 
that the method of uniqueness cannot be applied immediately. Instead one has to 
use integration by parts and other methods to produce vertices which satisfy the 
uniqueness condition before that integration technique can be used. Nevertheless 
the quark anomalous dimension and quark mass anomalous dimensions have been 
determined at 0(1/ Nj). For the quark anomalous dimension in arbitrary covariant 
gauge we have, 

[(2/x-l)(/x- 2) +a/j]C F r]o 



m 



[(2/i-l)(M-2)2V 



(5.10) 



and 



% = 



- /i(2// + pa - 5/x + 2){n - 1) [§ 2 (/x) + $(/i)] C A 



(8// - 92// + 270// - 30 V 
*(m)Ca 



124/1 - 12)((2|U - 1)(m- 2) + fjta) 



2(2/*- l)(2/z- 
3/t(/t — 1)\ploiCa 
[(32// 



3)(m-2; 

-2(2// + 



8//a 2 - 224//<S 



96// 

704// 5 - 3360// + 16//a 3 - 
6240// - 19//a 3 - 387// <S 2 
6// a 3 + 230/Za 2 - 222//a • 
Afj,a - 908// + 88)C A /i 
8(4// + 4//a - 32/i 4 - 13// a 



5/i 
6^ i 



2)C F ]9(/x) 



912// - 4//S 3 - 84//a 



,5^,2 



278//a 2 - 1124//a 
f 846// a - 6292// 
- 3416// - 48/ia 2 



+ 32/i 
/[4(2/x 



6)(2 M -l)(2/z-3)(/ i 
l)(2/ 1 -3)( M -l)(M-2)/i] 



+ 75// 
2)C F ] 



8//a - 70// - 2/ta 



2(2/ i -l/( At -2)2T 2 



where 



V 



Vi 



oo 
i=l f 



(5.11) 



(5.12) 



The expression for the quark mass anomalous dimension is independent of a and 
in the same notation, 

2Cf7/q 

^ - " (/, - 2)T> (5 6) 
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and, [55] , 

Lf~A, x | (4// 3 -13// 2 + 9/x~3) ^ 
^,2 = - |2(30(fO + ^-1)2 ) Cf 

(8// 5 - 92// 4 + 270m 3 ~ 301// 2 + 124// - 12)f (//) 

2//(/t - l)(2/i - 1)(2// - 3)(// - 2) 
- [16/t 6 - 128// 5 + 480// 4 - 900/i 3 + 831// 2 

- 344// + 44)]/[4//(2/i - 1)(2// - 3)(/z - 1) 2 (// - 2)] 

-^,)-t M )c] ^^:^ . ,5.4, 

Having focused on the operator associated with the coupling constant of QCD, 
one can also repeat the same analysis for the analogous operator of the non-abelian 
Thirring model which is the dimension two gluon mass operator. Whilst it clearly 
is not a gauge invariant operator its anomalous dimension has an interesting struc- 
ture. Inserting \A a ^A a ^ into the same Green's function as G®„G a ^, the 0{l/N f ) 
exponent in the Landau gauge is, |34j . 

™ = -mM^ + (^)- < 5 ' 15) 

Interestingly this can be rewritten as 

lA 2 (5c) = 7a(Sc) + lc{g c ) (5.16) 

in all dimensions d. It turns out that this is a general property of this operator in 
the Landau gauge. In j^l] it was shown that in QCD 

lA^(g) = ta(s) + 7c(ff) (5.17) 

to all orders in perturbation theory having first been shown at three loops in MS 
by explicit computation, |36|. The explicit four loop value is also now available in 
SU(N C ), |37|. Although this operator is gauge variant, it has been the subject of 
intense study in recent years as an effective gluon mass term. See, for example, 
|38| . Moreover, it has been studied in other gauges such as the maximal abelian 
gauge (MAG), |391 14UI ITS] . In this latter gauge the gauge field is written in terms 
of its diagonal (centre) and off-diagonal fields 

A^T A = A a ^T a + A^T* (5.18) 

where [T\T j ] = 0, T l £ {centre}, 1 < i < N%, 1 < a < NX and 1 < A < N A with TVf 
and iV| the respective dimensions of the centre of the group and its complement. 
(Here and in the next section we use the index A to denote the whole colour group, 
using the notation usually used in discussing the maximal abelian gauge.) Then in 
this gauge it turns out that the off-diagonal BRST invariant mass operator 

O = \AlA ai1 + ac a c a (5.19) 

satisfies a Slavnov- Taylor identity, similar to that for the analogous operator in the 
Landau gauge, which is, |39|. 

io{g) = lA-(g) -7c*G?) (5.20) 

where c* is the ghost in the centre of the colour group. Whilst this identity has 
been established on general grounds, for practical calculations, such as studying 
the condensation of the operator O of (|5.19|l . one needs the explicit values of the 
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anomalous dimensions. This requires the application of much of the earlier general 
discussion on renormalization and the maximal abelian gauge provides a compre- 
hensive example to illustrate these general remarks. 



6 QCD in non-linear gauges 

Given this interest in studying QCD in covariant gauges other than the canon- 
ical linear ones, we now focus on the explicit renormalization in the MAG. First, 
we define the gauge fixing for the maximal abelian gauge by recalling that for the 
Landau gauge. With l|5.18|l the gauge fixing is given by 



L Landau = §§ 



\a a A a » + \ac A c A 



(6.1) 



By contrast, the maximal abelian gauge is such that one minimizes the quadratic 
form in the off-diagonal sector only by fixing, . 



L^ AG = 88 



—A a A a ^ 
2 M 



s [?d»Ai] 



(6.2) 



where again 8 and 8 are the respective BRST and anti-BRST variations with 



sa; 

Sc a 
Sc a 

8c* 
Sc 1 



(V Q + !lf" J --\;y + 9f abc Alc c + gf abk Alc k ) 



1 



= gj abk c b c k + ^f abc c b c c 

= b a , 5b a = 

= - (d^ + gf^Ay) 

= \gf lhc c h c c 

= b* , 5b* = 



(6.3) 



and 



sa; 

Sc a 
Sc a 
Sb a 



- ( V° + 9 r c Ay + gr hc A\c c + gf abk Ay) 

- b a + gf abc c b c c + gf abk c b c k + gf abk c b c k 
1 



= 9f° 



'<■"<■" + t;<iJ '■"<■ 



g fabc b b-c _ g fabk b b-k + g fabk^b b k 

(dy + gf ibc Ay) , Sc l = -V+gf lbc c b c c 



^gf bc c b c c , 8b l = - gf lbc b b c c 



(6.4) 



The final term in l|6.2|l is required to avoid a singular ghost propagator. The residual 
gauge freedom in the diagonal sector is fixed by applying the usual Landau gauge. 
Though one can introduce an extra gauge parameter, a, for the inversion to obtain 
a centre gluon propagator. Varying 8 and 8 gives contributions to the interaction 
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Lagrangian with the remaining terms derived from (G A ) 2 . Thus, 



L gf - 



abk Ad -k 



cabk 



d»A a u c b c k 



1 

b f 
i 



+ c a d»d^c a + ffd^d^c 
1 



abc A a-b d n c c 



-s 



2/ 



d^A^AtA 



abk A k-a d t,-b 



k v 



cabk 



d"A k c b c c 



jacbd A a A bfj,-c c d jakbl ^ja j^b fi ^k j^l v _j_ jadcj j^a A j [i-c Q d 



2a 



_ jajcd A a A j MgC c d _|_ jajcl j^j MgC^ _|_ falcjj^aj^) ' PjffJ 



2 JO ' M 



i racbd-a -b c d 
+ g/o C C C C 



/, 

" fabcl 

A Jo 



4 

a 



abcd c a c b c c c d 



abcd c a c b c c c d 



-a-o c i i _pG 

c c c c + -/ 



c c c c 



f" 

•I O 



albc^a-b c c c l 



-f" 

i J o 



akbL c a c b c k c L 



9 



(6.5) 



where f^ 300 = f Am f cm and f ABCD = f ABe f CDe . Although this is a much 
more involved Lagrangian than that which results in the Landau gauge where there 
are a handful of interactions, it can be shown that it contains the usual covariant 
gauge fixed Lagrangian, 39 . Moreover, in the limit where the fields deriving from 
the centre of the colour group are set to zero one recovers the non-linear Curci- 
Ferrari gauge, . In other words the off-diagonal sector of the maximal abelian 
gauge corresponds to QCD fixed in the Curci-Ferrari gauge. Clearly the Lagrangian 
includes quartic ghost interactions which are always a feature of a non-linear gauge 
fixing. Though they do not destroy the renormalizability of the maximal abelian 
gauge Lagrangian which has been established by the algebraic renormalization tech- 
nology, [391 l4Hj . This is crucial to formulating the renormalization of QCD in the 
maximal abelian gauge as it produces the Slavnov- Taylor identities originating from 
the BRST transformations. More importantly it determines the form of the renor- 
malization of the fields and parameters as 



A ap. 



I ry a —a 




Za A a M , Ay* = ^Z^A 1 ^ 
c 1 



Z^ip , g a = H € Z g g 



(6.6) 



Consequently, one observes that the diagonal ghost 2-point function is finite. If this 
had not been deduced via algebraic renormalization, then in an explicit renormal- 
ization of the diagonal ghost 2-point function one would have erroneoulsy deduced 
that Z c i was unity. Instead to determine the correct Z c % one has to renormalize 
a Green's function containing one Z c i where all the other associated renormaliza- 
tion constants are known. For this case one choice would be the A^^c 1 vertex 
renormalization. Another useful feature is that the diagonal gluon wave function 
renormalization constant is not independent but is in fact related to that of the 
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coupling constant. Hence, the /3- function can be determined without resort to the 
renormalization of a vertex function. Moreover, this simplification is similar to what 
happens in QCD fixed in the background field gauge, [121 ■ Though in practical cal- 
culations it reduces substantially the number of diagrams to be determined as well 
as computation time which can be significant at high loop orders. The key point is 
that the preliminary analysis from algebraic renormalization not only determines 
the structure of the renormalization constants consistent with the underlying sym- 
metry, but also provides an efficient route for computing the anomalous dimensions 
themselves. 

Given that the maximal abelian gauge gluon mass operator does not renormal- 
ize independently, to have its explicit anomalous dimension requires only knowledge 
of the diagonal ghost anomalous dimension. Clearly given the large number of in- 
teractions in H6.5J1 , a three loop renormalization can only proceed with an automatic 
computer algebra approach for which the Mincer algorithm is the ideal tool, us- 
ing the Form version, and Qgraf, [21]. In completing the full three loop 
maximal abelian gauge renormalization of QCD, 37322 Feynman diagrams had to 
be considered, ^HJ. Though not all of these are non-zero. Some vanish trivially by 
the group theory structure of a diagram. For instance, where one has a self-energy 
insertion in a propagator line with a diagonal and off-diagonal field as external to 
this subgraph, then it will vanish due to 

fjk = fijc _ o ^ pbc ^ q ^ jabc ^ q _ ? ) 

The result of the renormalization procedure yields the anomalous dimensions, |18| . 
which have the structure illustrated by the example, 

7o(o) = ^[NZ((-3a + 35)C A -16T f N f ) + Ni((-6a-18)C A )]a 



' Nf ((- 6a 2 - 66a + 898)Ci - 560C A T f N f 



96NX 2 



384C F T f N f ) + N%N% ((- 54a 2 - 354a - 323)^ + 160C A T f N f ) 



+ 



n a ((- 60 " 2 - 372a + 510)C^) 
1 

6912iV| 3 



« 2 



Nf{{- 162a 3 - 2727a 2 - 2592C(3)a - 18036a 
- 1944C(3) + 302428)Ci + (6912a + 62208C(3) - 356032) C%T F N f 



+ (- 82944C(3) - 79680)C A C F T F N f + 4Q4Q8C A T F Nf 
+ 13824C F T F N f + 33792C F T F N?) 

+ NfN%((- 2754a 3 + 648a 2 C(3) - 28917a 2 - 4212aC(3) - 69309a 

+ 37260C(3) - 64544)Ci + (25488a + 103680C(3) 

- U072)C 2 A T F N f + (- 165888C(3) + 155520)C A C F T F N f 

+ \792QC A T F Nf) + N£Nf((- 7884a 3 + 22680a 2 C(3) - 84564a 2 
+ 97524aC(3) - 47142a + 433836C(3) - 56430)C^ 
+ (25056a - 124416<C(3) - \8UA)C A T F N f ) 

+ Nf((- 6480a 3 + 34992a 2 C(3) - 70092a 2 + 8424aC(3) 

+ 114912a +77112C(3)- 161028)Cl)] a 3 + (9(a 4 ) (6.8) 
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where a = g 2 /(167r 2 ) is conventionally used as the coupling constant. To be con- 
fident that the anomalous dimensions are in fact correct, the standard internal 
renormalization group checks we discussed above have been shown to hold. More- 
over, since the /3-function is a gauge independent object then the cancellation of a 
and the parameters and ATJ, provided not only a strong check on the calculation 
but on the group theory Form module which handles the consequences of splitting 
the gauge group into two sectors. Additional checks for this specific computation 
are provided by noting that the anomalous dimensions of the off-diagonal fields are 
correctly equivalent for all a to those of QCD fixed in the non-linear Curci-Ferrari 
gauge, which was introduced in 0T], in the limit N^/N^ — > where the explicit 
results were given in 03] EH3 ■ 



7 Renormalization and non-locality 

We now turn to more recent aspects of renormalization and that is the renor- 
malization of quantum field theories where a non- locality is present. The general 
properties of renormalization theory as discussed, for example, in 01], is based on 
the assumption that the Lagrangian is local. In other words there are no non-local 
interactions or operators. However, in certain field theories of interest, such as 
QCD, one encounters important operators which are non-local but whose proper- 
ties require investigation. We now summarize the status of two such studies where 
although there is a degree of non- locality present, it falls into a class which does 
not prevent calculations from being performed. By this we mean that the non-local 
operators can be rewritten in terms of a finite number of local fields and opera- 
tors. We refer to this as a localizable non-locality. By contrast, there are operators 
which, whilst being non-local, do not allow for a finite number of auxiliary fields to 
lead to a local Lagrangian. We term this a non-localizable non-locality. 

The first example is that of the Gribov problem which relates to the impos- 
sibility of globally fixing the gauge in a non-abelian gauge theory. This was first 
pointed out by Gribov, 0H], who proceeded to construct a path integral to study 
the problem. There the non-abelian gauge field is restricted to the region defined 
by the first zero of the Faddeev-Popov operator. This region, known as the Gribov 
volume, is of finite size and characterised by the Gribov mass, 7. The resulting 
path integral, [15]. 

Z = J DA 8{d»Al) Act {-d u D a v ) e~ s (7.1) 
essentially leads to a non-local Lagrangian 

L = — -G a ay G a ^ + Qlt A a_J_ A a lt _ (72) 

4 ^ 2 d v D u 2g 2 



The parameter 7 is not independent and satisfies the Gribov gap equation which 
is, at one loop, |45j. 



1 = C A 



3, (C A1 



-In 



V M 4 



0(a 2 ) . (7.3) 



Subsequently, in 0|)] BZ| Zwanziger managed to localize the original Gribov path 
integral by introducing a set of extra ghost fields {</>jf, <^ b , oj® b , u>^ b } where the last 
two fields are anti-commuting. This led to the Gribov-Zwanziger Lagrangian in the 
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Landau gauge, (IBIHT) . 

L GZ = L QCB + frbfigv (D v <j )fJ ) ab - u! ab ^d v {D v U), 1 ) ah 
- gf abc d u ujl e (A,c)V c " 

_ -Xl ffubcjpujfc + fabc A a»jfic\ + , ? q 

which clearly involves only a finite number of local interactions. Not only is the 
Lagrangian local, it is also renormalizable, 0H1 149) . which implies that one can 
perform loop computations and study the implications the Gribov parameter has on 
the infrared structure of QCD. The renormalization structure of (|7.4|l is interesting 
in that in the Landau gauge the renormalization constants of the extra ghost fields 
are not independent with, |4"%ll4l?] . 

Z<j, — — Z c . (7.5) 

Moreover, the renormalization of the Gribov parameter is not independent in the 
Landau gauge, satisfying, [4*§| . 

7 7 (a) = ~ [7 A (a) + 7 c (a)] (7.6) 

which is similar to the renormalization of ^A^A afJ -. Moreover, the quark, gluon, 
Faddeev-Popov ghost anomalous dimensions and the /3-function are unaltered by 
the presence of 7 and the extra Zwanziger ghost fields in MS. In a general linear 
covariant gauge one has the additional two loop MS results for the anomalous 
dimensions that 

70(a) = lw{a) = 7c (a) 

7 7 (a) = (16T F N f -(35 + 3a)C A )-^ 

a 2 

+ (X92C F T F N f +280C A T F N f - (449 - 3a))— + C>(a 3 ) (7.7) 

192 

which was deduced using Mincer. 

Equipped with these properties one can compute the corrections to the one 
loop mass gap equation which Gribov originally derived. This corresponds to the 
horizon condition definition of 47 , which is equivalent to 

r bc {A^(x)<p b ;(x)) = ^0- (7.8) 

which can be evaluated order by order in perturbation theory using the propagators, 

[SlEZlEnj, 

(A^AU-p)) = [{p2 y b f CAii] P,Ap) 

fabc 2 

{A a Jp)(ji c {-p)) = - 1 i>(p) 

xaczbd fabe £cde^4 

(<tf<P)4?(-p)) = —r,»» + ^i + c^f ^ (7 ' 9) 
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where the gluon propagator is suppressed in the infrared. Hence, at two loops in 

MS, ESDI, 



1 = C A 



5_3 (Crf 



V M 4 



/2017 11097 



A V 768 



2048 256 s ^ ' 48 \ /i 4 



128 \ X \-V 



1137 V5C(2) 2 ° 5 - 2 



CaTfN, 



2560 
7 



512 



24 ^ V ; 12 V A 4 



1 A /C^ 7 4 
8 r I"?" 



2 ^ 

TT 



0(a 3 



(7.10) 



where s 2 = (2^3/9) Cl 2 (27r/3). The one and two loop corrections are evaluated 
using the vacuum bubble approach discussed earlier. Though we note that the 
basic two loop multiscale vacuum bubble integral /(mf , to 2 , to 2 ;), where 



Til 2 2\ 

7(TO 1 ,m 2 ,TO 3 ) 



1 



/,■/ 



[k 2 — to 2 ][Z 2 — to 2 .]^ — I) 2 — to 



(7.11) 



has to be determined for to 2 € {0, iVUaj 2 



/ C J 47 2 }- For instance, the quanti- 



ties S2 and -\/5C(2) arise from the finite parts of 7(«v / Ca7 2 , 2VCU7 2 , «VCyi7 2 ) and 
I(iv / C J 47 2 , iy/CAj 2 , —iy/CAj 2 ) respectively. 

One interesting consequence of this gap equation is that it ensures the Kugo- 
Ojima confinement condition, |51j . of Faddeev-Popov ghost propagator enhance- 
ment is preserved at two loops. If one considers the full ghost propagator to have 
the form 

p l [l + u(P )J 

where u{p 2 ) represents the radiative corrections, then the Kugo-Ojima confinement 
criterion is that u(0) = — 1, Hence, the ghost propagator behaves as l/(p 2 ) 2 
as p 2 — > 0. As w(p 2 ) corresponds to the loop corrections of the ghost 2-point func- 
tion, then computing u(p 2 ) in the vacuum bubble expansion to two loops one can 
check if the Kugo-Ojima criterion holds at this order. Using the Gribov-Zwanzigcr 
Lagrangian, l|7.4[) . it transpires that the gap equation emerges as the u(0) + 1 term 
at 0(p 2 ) so that u(0) = — 1 precisely at two loops, [50] . 

Finally, we discuss a more recent study of non-locality in QCD, |52]. Earlier 
we considered the gauge variant dimension two operator ±A® A afl given the current 
interest in it as a potential mass operator for the gluon. However, if one wishes to 
have a gauge invariant Lagrangian with gluon mass then this operator is excluded. 
Instead to preserve gauge invariance and insist on a mass operator, one has to allow 
for a non-local mass operator. In [22], it was pointed out that aside from the one 
usually associated with the (non-renormalizable) Stueckelberg term, 
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(7.13) 
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there is another independent operator which is 

G%-^G a ^ (7.14) 
which can be added to the gauge fixed Lagrangian as 

L = i gf - ^G; v ^G a ^ . (7.15) 

This is a localizable operator in the sense we defined earlier. Consequently, | 52) . 
it can be localized by introducing the ghost fields {-B" , 5" H", S"} where the 
last two are anti-commuting, to give 

L = L g{ + ^(B; v - B%)G a ^ + - A Bl v {D"D a B^) a - \h^{D° D a H^) a . 

. (7 : 16) 

This Lagrangian has been analysed by algebraic renormalization and whilst it is 
renormalizable, it is not multiplicativcly rcnormalizable since new quartic ghost 
terms are generated through quantum corrections, [52] • Moreover, the localized 
dimension three operator (B" — B®)G aflv operator itself mixed into the lower 
dimensional mass operators [B^B**" - H^ v H a ^ v ) and (B* v - B^) 2 . Given 
that l|7.16[) is renormalizable, the one loop anomalous dimensions have been deter- 
mined. Clearly the addition of such an operator ought not to affect the established 
anomalous dimensions of the gluon, Faddeev-Popov ghost and quark fields and ex- 
plicit calculations verify this. Essentially within the respective 2-point functions 
the contributions from the extra ghost fields cancel. Moreover, 

7fl(o) = 1h{o) = {a-3)C A a + 0(a 2 ) (7.17) 

and the anomalous dimension of the mass operator itself can be deduced by inserting 
the gauge invariant equivalent operator (B^ u — B°; iy )G af * L ' into the S° -gluon 2- 
point function. The explicit Mincer calculation yields the result, |52| . 

lB G (a) = - (yC74 - ^TpN^j a 2 + 0(a 3 ) (7.18) 

which is not only independent of the gauge parameter as it ought to be, but is equiv- 
alent to the one loop QCD /3-function. This is the same as the one loop anomalous 
dimension of higher dimensional operators with similar Lorentz structure, |53j . In 
other words the operators G° v G a > iU , D fl G^ a D p G al " T , D v G a a p D^ D v G a ° p and 
D fJj D„D <J GpgD fl D 1 ' D a G a pB all have same one loop anomalous dimension as the 
non-local operator. We complete this section by remarking that it would be inter- 
esting to determine the higher loop corrections to l|7.18|l . 

8 Conclusions 

We close with general remarks. First, in the initial part of the article we have 
reviewed the techniques and general structure associated with the renormalization 
of quantum field theories in the context of O(N) (j) 4 theory and QCD. One of the 
interesting aspects of renormalization is the rich interplay between the abstract and 
mundane exercise of evaluating renormalization group functions and the implica- 
tions such results have on the underlying physics. For example, the anomalous di- 
mensions when evaluated at a critical point of the /^-function lead to predictions for 
the scaling behaviour of Green's functions which can be measured in experiments. 
However, from a technical point of view the critical exponents, when studied in the 
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large N expansion, also actually complement the checking of the explicit perturba- 
tive results. For large loop order computations this plays a role in establishing the 
correctness of the result as well as providing new information on the as yet unde- 
termined terms in the series at several orders in large N. This can be important 
due to the fact that such high order computations can presently only be under- 
taken with intense use of computer algebra and symbolic manipulation techniques. 
The latter part of the article dealt with the application of rcnormalization to more 
current computations including the renormalization of QCD in non-linear gauges. 
Although this follows the application of established techniques, such calculations do 
provide additional checks on the already determined three loop gauge independent 
renormalization group functions such as the /3-function. Whilst much of the theory 
of renormalization is based on the assumption of locality of the initial Lagrangian, 
we have also touched on very recent calculations of operators which are non-local 
in structure. Although these fall in the class of localizable non-local operators, it is 
possible that such studies might open the possibility of studying problems in QCD 
where such operators are important in probing the infrared structure of the theory. 
This is more evident, for example, in the Gribov problem where the structure of 
the ghost propagators satisfies the Kugo-Ojima confinement condition at two loops. 
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